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Abstract. We find a possibility of a weak universality of spin-glass phase transitions 
in three-dimensional ± J models. The Ising, the XY and the Heisenberg models seem 
to undergo finite-temperature phase transitions with a ratio of the critical exponents 
7/1/ ~ 2.4. Evaluated critical exponents may explain corresponding experimental 
results. The analyses are based upon nonequilibrium relaxation from a paramagnetic 
state and finite-time scaling. 



1. Introduction 

A spin-glass (SG) phenomenon has been attracting great interest both theoretically and 
experimentally [T|. Applications now cover a wide range of interdisciplinary fields of 
statistical physics and informational physics, as treated in this special issue. However, 
many subjects are not well understood, in spite of efforts made over almost thirty years. 
One of these subjects is whether or not the SG transition of real materials can be 
explained by a simple random-bond spin model. 

Spins of many SG materials are well-approximated by the Heisenberg spins. The 
simplest theoretical model is the Heisenberg model with random nearest-neighbour 
interactions. However, numerical studies suggest that there is no finite-temperature 
SG transition in this model [3 El- Kawamura |H E] proposed the chirality mechanism 
in order to solve this discrepancy. The chiral-glass (CG) transition occurs without the 
SG order. A small but finite random anisotropy in the real materials mixes the chirality 
degrees of freedom and the spin degrees of freedom. This anisotropy effect induces the 
SG transition observed in the real materials. The scenario is based upon results that the 
SG transition does not occur in the isotropic model. However, Matsubara et al. jEllZlE] 
recalculated the domain-wall excess energy and the SG susceptibility, from which they 
suggested that the finite-temperature SG transition does possibly occur. Methods are 
quite similar to the previous ones [21 E|- Subtle differences in the analyses of the obtained 
data drew an opposite conclusion. 

f Present address: BIRDS SYSTEMS RESEARCH INSTITUTE INC., 1-6-15 Hirakawa-cho, Chiyoda- 
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The spin-glass problem is one of the most difficult subjects in computational 
physics. It can be a tough bench-mark test for a new numerical method. It may 
be applied to other complex systems, if successful in the spin-glass investigations. The 
difficulty is caused by serious slow dynamics. It requires many Monte Carlo steps to 
reach the equilibrium states. An observed quantity at each step has a strong correlation 
even after the equilibration. The system sizes which can be treated in the simulations are 
accordingly limited to very small ones, e.g., mostly a linear size is twenty or less in three 
dimensions. Size effects are generally stronger in the continuous spin systems because 
the spins are soft and the boundary effect propagates faster. Frustration and randomness 
also yield a considerable size effect. The system sizes treated previously in the studies of 
the Heisenberg SG models are too small to extrapolate to the thermodynamic limit. This 
is our motivation for reexamining the SG transition using the nonequilibrium relaxation 
(NER) method P EH HH 021 ECU d • 
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Figure 1. A schematic diagram to approach the thermodynamic limit. 



The difficulty mentioned above can be overcome by using the NER method. This 
method takes an opposite approach to the thermodynamic limit. Figure ^ schematically 
shows a comparison between the conventional equilibrium simulational method and the 
NER method. In the conventional method one takes the infinite time limit first by 
achieving the equilibrium states in finite sizes. The thermodynamic limit is taken by 
the finite-size scaling analysis of the obtained data. In the NER method we take the 
infinite size limit first by dealing with a very large system within a finite time range 
before the finite-size effect appears. Then, the finite-time scaling analysis fFA. fTH] is 
performed to obtain the thermodynamic properties. The cost of a simulation is in the 
same order of L d+Z for both methods. However, a coefficient factor in the NER method 
is much smaller than that in the conventional method. An observation time length in 
the NER method is sufficient if we can observe a beginning of a final relaxation to the 
equilibrium states (equilibrium relaxation). On the other hand, it is necessary to wait 
until the end of the equilibrium relaxation in the conventional method. The latter time 
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scale is typically 10 ~ 10 2 times longer than the former one in the spin glass models. 
(For example, x sg of L — 17 in figure|U(a) or x sg at T = 0.56 in figured (a).) Therefore, 
the NER method has an advantage over the conventional method by this factor. We use 
the residual computational time to enlarge the system size and to increase statistical 
accuracy. 

By using the NER method we have made it clear that the SG transition occurs in the 
Heisenberg model at the same finite temperature as the CG transition occurs [T7j . The 
estimated critical exponent 7 is consistent with the corresponding experimental result 
[TH] . The chirality mechanism is not necessary to explain the spin-glass experiments 
since the chirality trivially freezes if the spin freezes. However, one may question the 
use and the validity of the NER method in the spin-glass phenomenon. Therefore, 
we have corroborated our method by studying the Ising SG model. Many numerical 
investigations [HH EUl EB 1221 I2HJ yield consistent results on the existence of the SG 
transition, the critical temperature and the critical exponents. They are also consistent 
with the corresponding experimental results |24j . The NER method yields consistent 
results for a small number of simulations as discussed in sectional 

In this procedure we have found a possibility of a weak universality: a critical 
exponent divided by u, for example 7/z/, is common among models in a weak universality 
class. A ratio of the critical exponents 7 / ' zv appearing in a finite-time scaling analysis 
is found to be consistent between the Heisenberg model and the Ising model. We 
have verified that a ratio 7/z/ is also consistent by evaluating the dynamic exponent z 
alone. The analysis is expanded to the XY SG model and the value is also found to be 
consistent. These findings are quite surprising. We must reconsider the role of the spin 
dimensions and the distribution of the randomness in the SG phase transition. 

This paper is organised as follows: In section |21 the model and the method are 
explained. Descriptions of the procedure of the NER method and the finite-time scaling 
are given. In section El the results on the Ising model, the Heisenberg model and the XY 
model are shown. Then, the possibility of a weak universality is discussed. Section |U is 
devoted to a summary. 

2. Model and Method - Nonequilibrium relaxation 

A model treated in this paper is the nearest-neighbour ± J random-bond model, 

H = J2 J ij S i ■ S r (!) 

M 

A linear size of a lattice is denoted by L. Skewed periodic boundary conditions are 
imposed, i.e., total numbers of spins iV = LxLx(L + l). An interaction takes two 
values of +J and — J with the same probability. The temperature T is scaled by J. 

Spins are updated by a single-spin-flip algorithm. The Metropolis (M) update is 
used in all models, whereas the heat-bath (H) update [2j is used in the Heisenberg model. 
Physical quantities observed in our simulations are the SG susceptibility Xsg, the CG 



Weak universality of spin-glass transition 



4 



susceptibility Xc g and the Binder parameter in regard to the spin-glass transition g sg . 
These quantities are calculated through the overlap between real replicas. 

First, we rewrite the thermal average by an arithmetic mean over thermally 
equilibrium ensembles labelled by a as 

1 m i m x,y,z 

m 



& ■ Sj) = - E sr ■ = -EE (2) 

a=l " l a=l M 

The bracket (• ■ ■) denotes the thermal average and m denotes a number of ensembles. 
The index /i stands for three components of spins: x, y and z. This expression is 
substituted into the definition of the SG susceptibility: 

m x,y,z 



N 



m- 



E E (<# f 



(3) 



where q?@ v = (l/N)J2i S^S^J is an overlap between the /i component of a spin i on an 
ensemble a: S^J and the z/ component of the spin on an ensemble (3: S^J . The bracket 
[• ■ -] c denotes the configurational average. 

Here, we introduce the following real replicas. Each real replica takes the same 
random bond configuration and the different paramagnetic initial spin state. They 
are updated in parallel with different random number sequences. This procedure 
corresponds to quenching from an infinite temperature. The thermal ensembles are 
realized by these real replicas which approach different equilibrium states. Therefore, 
we replace the thermal average by the average over these real replicas as equation (j2J. 
The indices a and (3 now represent real replicas. We do not take into consideration a 
constant term which arises from the overlap between the same replica a = (3 and use 
the following expressions in the simulations. 
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A number of replicas m controls the precision of the thermal average. It is better 
to take a large value. We prepare eight or nine replicas for each bond configuration 
in this paper. The scalar chirality is defined by three neighbouring spins as C, 
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forming a square plaquette in the direction from the i site. Constants in a definition 



of g sg are: A — 3, B — 1 for the Ising model, A 
9 for the Heisenberg model. 



6,5 = 4 for the XY model and 
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Figure 2. A schematic flow of our simulation. Solid bonds and broken bonds in the 
lattice depict ferromagnetic bonds and antifcrromagnetic bonds. 



Figure El shows a schematic diagram of the simulation procedure. We calculate 
a physical quantity at each time step t and obtain a relaxation function. Another 
simulation starts by changing a random bond configuration, initial spin states and a 
random number sequence. Then, another relaxation function is obtained. Finally, we 
take an average of data at each step over these different Monte Carlo runs. It should be 
noted that the average is over independent data. It guarantees an absence of systematic 
error due to correlations of the observed quantity, which we usually encounter in the 
conventional Monte Carlo time average. The obtained raw relaxation function is utilised 
by the following finite-time scaling analysis. 

The most important point in the NER method is to exclude the finite-size effect 
from the raw relaxation function. The method is based upon taking the infinite-size limit 
first. If a relaxation function includes a finite-size effect, it exhibits converging behaviour 
because every finite system has a definite equilibrium state. This behaviour misleads 
us into thinking that the temperature is in the paramagnetic phase even though it is 
the critical temperature. Therefore, the critical temperature is always underestimated if 
the size is insufficient. We check the size effect by changing the lattice sizes and always 
confirm a time range in which the size can be considered as infinity. 

The SG susceptibility is expected to diverge at the critical temperature (T sg ) as 
Xsg(t) ~ t 7//zi/ jini- We obtain T sg , 7 and zv by the finite-time scaling analysis on 
the relaxation functions of Xsg(^) i n the paramagnetic phase. ^7] Since the initial spin 
configuration is completely random, Xsg(^ = 0) ~ 1. We start a set of simulations at 
a temperature T that is obviously in the paramagnetic phase. The relaxation function 
X S g(t) at this temperature increases with t but soon converges to a finite value. As the 
temperature is lowered to approach the critical temperature, the relaxation function 
tends to show diverging behaviour. Since the temperature is still in the paramagnetic 
phase, the relaxation finally converges to a finite value after a correlation time t(T). 
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The spin-glass correlation increases with time and reaches the correlation length £(T) 
after this correlation time. Two quantities relate with each other by z as r(T) ~ £, Z {T). 
Therefore, the correlation time should diverge at T sg as 

T(T)~(T-T sg )-*\ (7) 

The correlation time can be estimated by scaling the raw relaxation function. We obtain 
y/zv and t(T) so that the scaled functions Xsg(t)t~' 7 ^ zu at all temperatures plotted 
against t/r(T) fall onto a single curve. Then, the critical temperature and the exponent 
zv are estimated by the least-squares fitting with the equation (jZJ). Since a ratio y/zv 
is already estimated by the scaling, 7 is obtained. 

The NER of the Binder parameter g sg (t) is calculated at the obtained T sg . Since 
quantity is related to the fourth-order cumulant, many bond samples are necessary to 
obtain meaningful data. The number of bond configurations to obtain the results in 
this paper is summarised in tabled The Binder parameter is expected to diverge at T sg 
as g sg (t) x L d ~ t d / z [Hj, by which z is independently obtained. Then, v is estimated 
from a value of zv obtained by the r-fitting explained above. All exponents are now 
estimated by the scaling relation. It is possible to compare the critical exponents with 
the experimental results. 

The last procedure of our method is to corroborate the results by observing the 
NER of Xs g at the obtained T sg . It should diverge as t 1 ^" with the same exponent 
obtained by finite-time scaling. If the exponents are inconsistent, the scaling analysis is 
misled by an insufficient time range or by the finite-size effect. In the Ising model we 
perform another check at T sg by observing the NER of the distribution function of the 
replica overlap, P(q,t). The finite-time scaling plot of P(q,t) should ride on a single 
scaling function with the same exponent obtained by finite-time scaling of Xsg- This is 
a direct interpretation of the finite-size scaling of P(q, L) J2D] by t oc L z . 

3. Results 

Numbers of bond configurations to obtain data at the critical temperature are 
summarised in table d The numbers at other temperatures are mostly in the same 
order. For each bond configuration, we prepared eight replicas for the XY and the 
Heisenberg model and nine replicas for the Ising model. 

3.1. Ising model 

Figure |3] shows an analysis to determine the critical temperature and the exponent. 
The simulation is performed just to check that our method gives results consistent with 
previous investigations [El 1201 EH 1221 EHj - Therefore, the system size is very small 
(L = 19) and the time range is very short. Finite-size effects are found to appear for 
t > 5000 by comparing with results of L = 29. Only data before this time are used in 
the scaling analysis. Figure El (b) is an example of finite-time scaling. A choice of y/zv 
is possible for y/zv = 0.38 ~ 0.40. A set of the correlation time at each temperature is 
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Table 1. Numbers of bond configurations to obtain data of \sg and g sg at T sg in this 
paper. Indices (M) and (H) in the Heisenberg model denote update algorithms: (M) 
for the Metropolis and (H) for the heat-bath. Arrows mean that the number is same 
as to the right. 
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255480 85480 18576 12626 
-> -> -> 7803 
43114 18316 7038 



1830 172 



10" 



60 9 

" 10 



10 J 



7=1.375 

7=1.400 

7=1.425 
7=1.450 
7=1.475 
T= 1.500 



(a): 




10 J 



10 1 



1600 bond samples 



10 z 



10- 



10' 



10 M 



10 



T= 1.375 - 

r=i.400 

r=1.425 
7 = 1.450 

r=1.475 

7=1.500, 

10 2 



(b): 



y/zv=0.39 




T-T. 



10 3 10 4 



tKT) 



sg 



Figure 3. (a) The NER of Xsg of the Ising model at high temperatures, (b) The 
finite-time scaling plot for a choice of ~ijzv — 0.39. We obtain r(T) and ^jzv so that 
this scaling plot is good. The scaling is also possible for ^/zv = 0.38 ~ 0.40. (c) The 
least-squares fitting of t(T) supposing t(T) cx \T — T sg \~ zu . 



estimated for each choice of this exponent. Then, the critical temperature is obtained 
as summarised in table El As the exponent increases, T sg decreases. We ignore a result 
of "jjzv = 0.400 which deviates a lot from the others. Our estimates are 

T sg = 1.17(4), j/zv = 0.3875(75), zv = 9.3(12), 7 = 3.7(5). (8) 
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Table 2. A list of the critical temperature and an exponent zv obtained by the 
finite-time scaling analysis in the Ising model. A ratio of exponents 7/ zv denotes the 
possible value in the finite-time scaling. The least-squares fitting errors are denoted 
by x 2 - 



7/ ' zv 


T 


zv 


7 


x 2 


0.400 


1.05 


12.63(18) 


5.1(1) 


1.40 


0.395 


1.13 


10.07(37) 


4.0(2) 


0.66 


0.390 


1.20 


8.46(28) 


3.3(1) 


0.54 


0.385 


1.22 


8.67( 9) 


3.3(0) 


1.49 


0.380 


1.21 


9.82( 9) 


3.7(0) 


1.84 



The results of the finite-time scaling analysis are checked by the raw NER data 
at the obtained critical temperature. Figure 0] (a) shows relaxation data of x sg and 
g sg x L d . The SG susceptibility diverges algebraically with an exponent ^jzv = 0.38 
that is consistent with the scaling result ^y/zv = 0.3875(75). The critical relaxation 
process begins around t ~ 100 and seems to continue to infinity. The Binder parameter 
also shows diverging behaviour with t d ^ z , from which we obtain the dynamic exponent 
z = 6.2(2). Then, an exponent v is estimated as v = 1.5(3). A ratio of the critical 
exponents 7/V = 2.4(1). The obtained results are consistent with previous numerical 
investigations [THl I2D1 El 1221 12H] and the corresponding experimental results [21] as 
summarised in table 01 Since the lattice size and the time range are insufficient, the final 
numerical results have large error bars. As discussed in the previous section the critical 
temperature may be underestimated by using a small lattice. We plan to estimate them 
with high accuracy by large-scale NER analyses. 

A time evolution of the distribution function of the overlap P(q, t) at T = T sg = 1.17 
is shown in figure EJ^b). The system size L = 17. It exhibits a single Gaussian form with 
a peak at q = before the size effect of x sg appears at t = 10 5 as shown in figure|Ha). As 
the time increases, the width of the distribution grows in accordance with the divergence 
of the spin-glass susceptibility. It is possible to scale P(q,t)/t y ^ 2zi/ plotted versus 
for various time steps from t — 10 to t — 10 4 (figure 0] (c)). The critical exponent j/zv 
is also consistent with the finite-time scaling of Xsg- The scaled data deviate a little for 
t = 10 because the time is just before the relaxation of Xsg reaches the critical relaxation 
region as shown in figure |U (a). The distribution changes its shape to having two peaks 
at ±g eq after the finite-size effect appears. The shape is flat at this crossover time. 

It is found that the NER function knows the critical phenomenon from its very 
early time steps: t — 10 ~ 100. The NER method is now clearly shown to be applicable 
to the spin glass phenomenon. 

3.2. Heisenberg model 

We apply the same analysis performed in the Ising model to the Heisenberg model. 
Finite-time scaling results have already been shown briefly in Ref. ^7] and the detailed 
analysis will be reported elsewhere. The system size is L = 59 and the time scale is 
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Figure 4. (a) The NER of x sg and g sg x L d of the Ising model at T sg = 1.17. Two 
lines, Vl zv with j/zv = 0.38 and t d / z with z = 6.2, are guides for eyes, (b) The NER 
of the distribution function P(q,t) at T sg for L = 17. The shape changes from the 
single-peaked to the double-peaked when the size effect of \sg appears at t — 10 5 . (c) 
The finite-time scaling plot of P(q,t). (d) A three-dimensional plot of P(q,t). 

70000 Monte Carlo steps. Typical numbers of bond configurations are same as given in 
tabled The finite-time scaling results are 

T sg = 0.20(2), i/zu = 0.39(5), zv = 4.8(10), 7 = 1.9(5). (9) 

These results are checked by the raw NER at T = 0.21 as shown in figure The SG 
susceptibility diverges algebraically with an exponent 'j/zv = 0.38, which is consistent 
with the result of the finite-time scaling. We performed simulations of both Metropolis 
update and heat-bath update. The Metropolis result denoted by (M) and the heat- 
bath result denoted by (H) exhibit the same critical behaviour, while the amplitudes 
are different by a factor of 3.5. NER behaviours are independent from the update 
algorithm. The consistency supports the criticality of the SG order at this temperature. 
It is noted that the critical divergence begins at a very early time: t ~ 100. 

The Binder parameter exhibits a critical divergence t d/z with z = 6.2(5). The value 
is consistent with that in the Ising model. Since zv is obtained by finite-time scaling, v 
is estimated as v = 0.8(2). A ratio of the critical exponent 7/1/ = 2.3(3). The results 
are compared with an experimental result ^H] m table El They are not inconsistent. 
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Figure 5. The NER of Xsg and g sg x L d of the Heisenberg model at T sg = 0.21. 
Lattice sizes are L — 39 for g sg x L d , L — 89 for x sg (H) and L — 59 for x sg (M). Indices 
(H) and (M) denote the heat-bath updated and the Metropolis updated, respectively. 
Lines, Vl zu with ^jzv = 0.38 and t d l z with z = 6.2, are guides for eyes. 

3.3. XY model 

It has been considered that there is no SG transition in this model [23l2n]- Only the CG 
transition with respect to the vector chirality is considered to occur [23123122] • However, 
a possibility of the SG transition has recently been identified by several investigations 
[3UJ EU EH! We applied the NER analysis on this model and our result supports the 
latter conclusion: the SG transition occurs. 

Our finite-time scaling analysis on the XY model is not yet conclusive in regard 
to whether the SG transition and the CG transition occur at the same temperature or 
not. A system size [L = 39), Monte Carlo steps (10 5 ) and a temperature range used in 
the scaling analysis are insufficient to extract a conclusion. However, as we increase the 
size and the steps, both critical temperatures seem to approach each other: T sg increase 
from low and T cg decrease from high. Therefore, we consider that both transitions 
occur simultaneously. Investigations are now being carried out and the details will be 
reported elsewhere. What has now been made clear is that both transitions occur in 
a temperature range of 0.4 < T < 0.46. In this paper we do not examine the issue of 
simultaneous transition but focus on the existence of the SG transition. 

Figure Efa) shows raw NER plots of Xsg an d Xcg near and above the critical 
temperature. There is no difference in Xsg between T = 0.43 and T = 0.46 within 
the present time steps. They exhibit a critical divergence with the same exponent and 
amplitude. The SG transition is considered to occur near T = 0.43. From the slope 
we obtain an exponent 'j/zu = 0.35. Note that this value is a little smaller than that 
of the Ising model and the Heisenberg model (j/zv ~ 0.38). The NER of the Binder 
parameter is shown in figure Efb). It exhibits a critical divergence with an exponent d/z 
with z = 6.8(5), which is also a little larger than that of the other models. However, we 
obtain a ratio of the critical exponents 7/V = 2.4(2), which is consistent with the other 
models. 
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Figure 6. (a) The NER of Xsg in the XY model above the critical temperature 
T sg ~ 0.43. Lattice size is L = 39. A line f/ zv with 1 / zv = 0.35 is a guide for eyes, 
(b) The NER of the Binder parameter multiplied by L d . A line t d / z with z = 6.8 is a 
guide for eyes. 



3.4- Weak universality 

The SG transition occurs in all models as shown in the preceding subsections. A ratio 
of the critical exponents j/u takes a common value around 2.4. Therefore, there is a 
possibility of weak universality among these transitions. Not only a value of 7/V but 
the NER functions themselves suggest that the transitions are qualitatively equivalent. 

Figure [7| (a) shows the NER functions of Xsg at the critical temperature for all the 
models treated in this paper. The data of the Heisenberg model with the Metropolis 
update are multiplied by a factor 3.5 in order to compare with a result of the Ising 
model and that of the Heisenberg model with the heat-bath update. These three NER 
functions are not distinguishable. If we take into account a correction-to-scaling term, 
the relaxation functions can be fitted from the first few steps (Bold lines in figured (a)) 
by an expression [23 : 



Here, exponents of the leading term are set 7/1/ = 2.356 and z = 6.2. The correction- 
to-scaling exponent w = 3. Coefficient constants are A = 7.6 and B = 0.7. The same 
expression also fits the NER function of the XY model, but with the dynamic exponent 
z = 6.8 and a constant A = 3.3. The NER functions of the Binder parameter are 
shown in figure [7| (b). If we multiply the result of the Ising model by a factor 2.2, it is 
indistinguishable from the curve of the Heisenberg model. 

4. Summary 

By applying the nonequilibrium relaxation method it has been made clear that the ± J 
models in three dimensions undergo finite-temperature spin-glass transitions. There is 
a possibility that these models belong to the same weak universality class with a ratio of 
the critical exponents 7/1/ ~ 2.4. We compare our results with other numerical results 
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Figure 7. (a) NER plots of the Xsg at T sg . Correction-to-scaling fittings are depicted 
by bold lines with 7/1/ = 2.356, w = 3, z = 6.2 for the Ising/Hcisenberg model and 
z = 6.8 for the XY model. NER functions except the XY model are indistinguishable. 
The data of the Heisenberg model with the Metropolis update are multiplied by 3.5. 
(b) NER plots of the Binder parameter multiplied by L d . The relaxation function of 
the Ising model multiplied by 2.2 coincides with that of the Heisenberg model with the 
heat-bath update. 



and the experimental results in table 0J They agree well within the numerical errors. 
Since the error bars are rather large at present, further efforts to improve precision are 
necessary in order to prove weak universality. 

The spin-glass transition of the Heisenberg model is found to be very similar to 
that of the Ising model. The relaxation functions of Xsg and g sg and values of a ratio of 
the exponents 7/1/ and the dynamic exponent z are consistent between the two models. 
If one considers that the spin-glass transition occurs in the Ising model, it may be 
thought that it occurs in the Heisenberg model in the same accuracy. Only the dynamic 
exponent of the XY model differs from the other models. Spin-glass transition and weak 
universality in models with Gaussian bond distributions is a problem to be checked in 
future work. The type of bond distributions may be important. 

The NER method has been shown to be particularly effective in the spin-glass 
study. Critical behaviour is observed from very early time steps even though it takes 
a very long time to achieve the equilibrium states. What is long is the nonequilibrium 
relaxation process after a short initial relaxation before the final equilibrium relaxation. 
This long process is discarded in conventional simulations, while it is utilised in the 
NER method. This is one reason why the NER method is advantageous in this system. 
Applications to various complex systems with slow dynamics are fruitful [33J HH] . 
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Table 3. Estimates of T sg , 7, v, 7/1/ and z in the ±J models in three dimensions. 





T 


7 


V 


j/v 


z 


Ising SG 












Present work 


1.17(4) 


3.6(6) 


1.5(3) 


2.4(1) 


6.2(2) 


Ref. |2I] 


1.11(4) 


4.0(8) 


1.7(4) 


2.35(5) 




Ref. 123 


1.195(15) 


2.95(30) 


1.35(10) 


2.225(25) 


5.65(15) 


Experiment 24 j 




4.0(3) 


~ 1.7 


~ 2.4 




Heisenberg SG 












Present work 


0.20(2) 


1.9(5) 


0.8(2) 


2.3(3) 


6.2(5) 


Ref. |H] 


0.18(1) 


2.0(2) 


0.97(5) 


2.1(1) 




Experiment [TSJ 




2.3(4) 


1.25(25) 


2.0(7) 




XY SG 
Present work 


0.43(3) 






2.4(2) 


6.8(5) 
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